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Abstract. The structure and energetics of superflow around quantized vortices, and 
the motion inherited by these vortices from this superflow, are explored in the general 
setting of the superfluidity of helium-four in arbitrary dimensions. The vortices may 
be idealized as objects of co-dimension two, such as one-dimensional loops and two- 
dimensional closed surfaces, respectively, in the cases of three- and four-dimensional 
superfluidity. By using the analogy between vorticial superflow and Ampere-Maxwell 
magnetostatics, the equilibrium superflow containing any specified collection of vortices 
is constructed. The energy of the superflow is found to take on a simple form for 
vortices that are smooth and asymptotically large, compared with the vortex core 
size. The motion of vortices is analyzed in general, as well as for the special cases of 
hyper-spherical and weakly distorted hyper-planar vortices. In all dimensions, vortex 
motion reflects vortex geometry. In dimension four and higher, this includes not only 
extrinsic but also intrinsic aspects of the vortex shape, which enter via the first and 
second fundamental forms of classical geometry. For hyper-spherical vortices, which 
generalize the vortex rings of three dimensional superfluidity, the energy-momentum 
relation is determined. Simple scaling arguments recover the essential features of these 
results, up to numerical and logarithmic factors. 



PACS numbers: 67.25.dk,67.25.dg 
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1. Introduction and aims 

First conceived of and analyzed by Onsager [H [3] and Feynman [H [5] , quantized vortices 
are topologically stable excitations of superfluid helium-four. They are manifestations of 
quantum mechanics which owe their stability to the quantization law that the circulation 
of superflow obeys, leaving their imprint at arbitrarily large distances. Vortex excitations 
play an essential role in determining the equilibrium properties of superfluid helium-four, 
especially in the mechanism via which superfluidity is lost as the temperature is raised 
through the "lambda" transition to the normal fluid state. They also play an essential 
role in the mechanism via which superflow in toroid-shaped samples is dissipated. For 
an in-depth account of the properties of quantized vortices, see Ref. [6j. 

When the space filled by the superfluid helium-four (i.e., the ambient space) is two- 
or three-dimensional, and when they are viewed at length-scales larger than their core 
size, quantized vortices are, respectively, the familiar, geometrical point-like vortices 
of two-dimensional superfluidity or line-like vortices of three-dimensional superfluidity. 
The aim of the present paper is to examine some of the characteristic properties of 
quantized vortices that are familiar in two- and three-dimensional superfluid helium-four, 
and to determine how these topological excitations and their characteristic properties 
can be extended to situations in which the superfluid fills spaces having arbitrary 
numbers of dimensions. For the case of three dimensions, these properties are analyzed in 
Ref. [7] . Properties that we shall address here include the structure of quantized vortices 
and the flows around them, the motion that vortices inherit from such flows, and the 
flow energetics, all in a variety of settings. In all dimensions, these excitations (again, 
when viewed at length-scales larger than their core size) continue to be geometrical 
structures having codimension two: their dimensionality grows with the dimension of 
the ambient space, always remaining two dimensions behind. Nevertheless, we shall use 
the terminology "quantized vortex" (or, for simplicity, just "vortex") for such excitations 
in all dimensions, not only two and three. 

We have chosen to present our developments using the traditional language of 
vector and tensor calculus. However, along the way we shall pause to mention how 
these developments appear when formulated in the language of exterior calculus and 
differential forms; see Ref. [H]. This language is in fact tailor-made for the questions at 
hand, and enables a highly economical development. However, we have not relied on it, 
as there may be readers for whom it would require an unfamiliar level of abstraction. 

This paper is organized as follows. In Section [21 we discuss the structure of vortices 
in superfluid helium-four in arbitrary dimensions of space, first from the perspectives 
of length-scales and topology, and then by developing the conditions that flows of a 
specific vorticial content must obey, in both the integral and the differential versions. 
In Section [3l we introduce considerations of kinetic energy, examining superflows of 
given vorticial content at equilibrium, and solving for the associated velocity fields 
of these flows. In Section HI we analyze the dynamics that vortices inherit from the 
equilibrium flows associated with them. Here, we specialize to asymptotically large. 
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smooth vortices, for which the dynamics turns out to be a simple and clear reflection 
of the intrinsic and extrinsic geometry of the shapes of the vortices. To illustrate our 
results, we consider the velocities of spherical and hyper-spherical vortices, as well as 
the small-amplitude excitations of planar and hyper-planar vortices. In Section [5l 
we return to issues of energetics, and determine the equilibrium kinetic energy of 
superfiows containing individual vortices in terms of the geometry of the shapes of the 
vortices. We illustrate this result for large, smooth — but otherwise arbitrarily shaped — 
vortices, for which a simple, asymptotic formula holds, and we give an elementary 
argument for this asymptotic formula. Here, we also address the energetics of maximally 
symmetrically shaped spherical and hyper-spherical vortices, and we use our knowledge 
of the energies and velocities of the spherical and hyper-spherical vortices to determine 
their momenta and, hence, their energy-momentum relations. In Section [6], we use 
scaling and dimensional analysis to obtain the aforementioned results, up to numerical 
factors and logarithmic dependencies on the short- distance cut-off, i.e., the core size 
of the vortices. We end, in Section U\ with some concluding remarks. We have 
streamlined the presentation by relegating to a sequence of appendices much of the 
technical material. 

2. Structure of superfiuid vortices in arbitrary dimensions 

2.1. General morphology of vortices 

We consider a large, D-dimensional, Euclidian, hyper-cubic volume Q of superfiuid, with 
each point in Q being specified by its D-component position-vector x, the Cartesian 
components of which we denote by {xci}d=i- the faces of fl we suppose that all fields 
obey periodic boundary conditions. We shall be concerned with flows of superfiuid in Q, 
and these we describe via their D- vector flow fields V{x). The phenomena that we shall 
be addressing occur on length-scales larger than the characteristic superfiuid coherence- 
or healing- length ^, which determines the length-scale over which the superfiuid density 
falls to zero as the centre of a vortex is approached (i.e., the vortex core size). 

Provided we do indeed examine flows on length-scales larger than ^, and hence 
restrict our attention to long-wavelength, low-energy excitations of the superfiuid, the 
amplitude {"^{x)] of the complex scalar superfiuid order parameter field \E'(x) can be 
taken to be a nonzero constant |\&o| away from the vortices. In addition, these vortices 
can be idealized as being supported on vanishingly thin sub-manifolds of the ambient 
D-dimensional space. The remaining freedom lies in the phase factor \E'(x)/|\E'(a;)| of 
"^{x), and — except on vortex sub-manifolds, where it is not defined — this takes values 
on the unit circle in the complex plane. Owing to the topology of the coset space 
(in this case, the unit circle) in which this phase factor resides, and specifically the 
fact that its fundamental homotopy group Hi is the additive group Z of the integers 
n = 0,±1,..., these vortices are indeed topologically stable defects in the superfiuid 
order [9j. The strength n of any vortex is the number of times the phase factor winds 
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through 271 as the vortex is encircled a single time. For D-dimensional superfluids, and 
at the length-scales of interest, these vortices are supported on manifolds of dimension 
D — 2, a. combination that will occur frequently, so we shall denote it by Ip. Said another 
way, the vortices reside on co- dimension- 2 sub-manifolds of Z)-dimensional space, which 
means that we would have to add two further dimensions to "fatten" each vortex into 
a space-filling object. In any dimension D > 2, there are one- dimensional closed paths 
that loop around the vortices. As we shall discuss shortly, the circulation of the flow 
around such loops is quantized and nonzero, whereas the circulation around closed paths 
that do not loop around a vortex is zero. 



In the familiar case of superfluids in three dimensions, the vortices are geometrical 
objects that are supported on one-dimensional sub-manifolds (see Fig. [1]). The loops 
that encircle the vortices, along which the circulation of the flow is quantized, are also 
one-dimensional objects In the equally familiar case of superfluids in D = 2, each 
vortex is a pair of opposing-circulation zero- dimensional objects (i.e., a structure that 
it usually called a vortex/anti- vortex pair). The loops around one or other member 
of the pair, however, along which the circulation of the flow is quantized, remain one 
dimensional. (Strictly speaking, in D = 2 each vortex is supported on a pair of points, 
which is not a manifold, but we shall admit this abuse of nomenclature.) 

Our main objectives are to explore the spatial structure and energetics of flows 
possessing a set of vortices at prescribed locations, as well as the motion that such 
flows confer on the vortices themselves. The flows that we shall consider are states of 
thermodynamic equilibrium, but are maintained away from the no-flow state, V = 0, 
because they are constrained to have a prescribed vorticial content. 

I The vortices themselves are closed loops that may or may not wind through the boundaries of the 
sample, a possibility allowed by the Z?-torus topology of the sample. This topology arises from the 
adopted periodic boundary conditions but will not play a significant role here. 




Figure 1. Vortices in three-dimensional superfluids. They form closed one- 
dimensional loops, possibly via the periodic boundary conditions. 
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2.2. Fixing the vorticial content of superflows 

The circulation n of the flow along any one- dimensional closed contour is defined to be 
the line-integral of the corresponding velocity field along that contour: 



K 



Vd{x)dxd. (1) 



Here and throughout this paper, a summation from 1 to D is implied over repeated 
indices such as d, di, etc. The vorticial content of the flow can then be prescribed 
via the following condition. For a vortex of strength n, the circulation k, is quantized 
according to the number of times TV the contour encircles the [i.e., {D — 2)]-dimen- 
sional sub-manifold supporting the vortex: 

K = {2nh/M)Afn, (2) 

where h is Planck's constant and M is the mass of one of the particles whose 
condensation causes the superfluidity. We adopt units in which h/M has the value 
unity. In the language of differential forms, the flow field is described by the 1-form 
V{x), and the quantization of circulation ([2]) is expressed as 

V = 2n^^n. (3) 

We shall be concerned with vortices of unit strength (i.e., n = 1). Owing to the linearity 
of the theory, results for higher-strength vortices can be straightforwardly obtained 
from those for unit-strength vortices via elementary scalings with n (i.e., velocities scale 
linearly with n; energies scale quadratically) . 

For a flow V to have the prescribed vorticial content, it must be expressible as a 
sum of two contributions: 

(i) a real vector field V{x) that is singular on any vortex-supporting sub- manifolds, and 
has the appropriate quantized circulation around them, and thus its contribution 
to the flow includes the flow's vorticial part; and 

(ii) the gradient of a real, scalar phase field $(x) that is smooth and single-valued 
throughout Q. Owing to the latter property, this component of the flow is capable 
of contributing only to the irrotational (i.e., circulation-free) aspects of the flow. 

Thus, we consider flows of the form 

Vd{x) = Vd{x) + VdHx), (4) 

in which the gradient denotes the partial derivative d/dxd with respect to the d^^ 
spatial coordinate Xd- (The Vd$ term carries a dimensional factor of h/M, which we 
set to unity, earlier.) In the language of differential forms, Eq. (jlj) is expressed as 

V = y + d$, (5) 

where d^ is the exact 1-form that results from taking the exterior derivative of the 
0-form (i.e., the function) $. 
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To ensure that the vorticial part V of the flow V has the appropriate circulation 
on loops surrounding any vortex-supporting sub-manifolds, it must obey the following 
inhomogeneous partial differential equation: 

ed,-dn^do-iVdnix) =2nJd,...d^{x), (6) 

where ed^...cip is the completely skew-symmetric D-dimensional Levi-Civita symbol. 
Equation ([6]) is a natural generalization to D-dimensions of the static version of 
the three-dimensional differential expression of the Ampere-Maxwell law relating the 
magnetic field (the analogue of V) to the electric current density (the analogue of J), 
which in suitable units reads 

curlV = 27rJ. (7) 

An account of the analogy between superfluidity and magnetostatics for the case of 
three dimensions is given in Ref. [TO] . 

The source field in Eq. ([6]), J, is a completely skew-symmetric, rank ij) tensor field, 
which obeys the divergencelessness condition dd^Jd^...dj^{x) = 0, i.e., the consistency 
condition that follows, e.g., from the application of dd^ to Eq. ([6]). (We discuss this 



point further in Appendix A ) This field J encodes the singularity in the density of 
the vorticity associated with any ^-dimensional vortices in the superflow. The factor 
of 27T in Eq. ([6]) is extracted to ensure that unit-strength sources gives rise to unit 
quanta of circulation. That V must obey Eq. ([6]) follows from the integral form of 
this law, Eq. ([21), which can be obtained from Eq. ([6]) via the D-dimensional version of 
Stokes' theorem. The right-hand side of Eq. ([2]) is the flux of the source J{x) through 
the corresponding closed contour. (In its integral form, the three-dimensional Ampere- 
Maxwell law dictates that the line-integral of the magnetic field along a loop in three- 
dimensional space has a value determined by the flux of the electric current threading 
the loop.) The skew-symmetric tensor nature of the source field J, present for D > 4, 
is a natural generalization of the vectorial nature of J present in D = 3, and reflects the 
corresponding higher-than-one-dimensional geometry of the sub-manifold that supports 
the vortex, as we shall discuss shortly. 

In the language of differential forms, Eq. ([6]) is expressed as an equality of ^-forms, 

i.e., 

dV = 271 J, (8) 

where dV is the 2-form that results from taking the exterior derivative of V, the form 
-k dV is the ^-form dual to it that results from applying the (Eucidian) Hodge-T*r operator 
to it, and J is the source ^-form. The dual of this equation, which is an equality of 
2- forms obtained by applying T*r to it, reads 

dV = ^{27Tj). (9) 

In the language of differential forms, the divergencelessness of the source J, mentioned 
shortly after Eq. amounts to J being a closed form, i.e., dJ = 0. 

What value should the source field J have if the D-dimensional flow is to contain 
a single, unit-strength, quantized vortex? Recall that for D = 3 (so that Ip = 1), the 
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vortex is one-dimensional and supported on the curve x = R(si), parameterized by the 
single variable Si (see Fig. [2]). In this case, the source takes the value 



associated with the analog of a unit-strength electrical current flowing along the curve 
and supported on it. The sign of this analogue of the current (i.e., its direction around 
the loop) determines the sense in which the superflow swirls around the vortex core, 
via a right-hand rule. Observe that the geometry of the line enters not only in terms 



of the curve R(si) on which the vortex is supported, but also via the tangent vector 
to this curve, dH/dsi. The geometrically natural extension of the source to arbitrary 
dimension D involves: 

(i) concentrating the source on the appropriate ^-dimensional sub-manifold; 

(ii) specifying the shape of this sub- manifold, via x = R{s), where s denotes the set of 
the Ip variables {si,S2, . . . ,Sjp} that, when they explore their domains, trace out 
the idealized core of the vortex; and 

(iii) endowing the source with appropriate geometrical content, via the natural 
generalization of the tangent vector. 

This generalization is provided by the oriented, ^-dimensional (tensorial) volume 
^ai-a^idaj^Rdi) ■ ■ ■ (daj^Rd^) (whcrc da dcuotes the partial derivative with respect to 
the a^^ independent variable that parameterizes the vortex sub-manifold, in this case, 
d/dsa) of the parallelepiped constructed from {dR/ dsa}a=ii i-^-; ^ linearly 

independent, coordinate tangent- vectors to the sub-manifold at R{s). Together with 
the measure d^s, this volume plays the role of an inflnitesimal segment of vortex line 
(in D = 3) or patch of vortex area (in D = 4, as sketched in Fig. [3]). Here, ^ai-. a^ is 
the completely skew-symmetric ^-dimensional Levi-Civita symbol, in contrast with the 
D-dimensional 6^^...^^, which we deflned shortly after Eq. ([6]). 




(10) 




Figure 2. A loop vortex in a three-dimensional superfluid specified by R{s) with 
arc- length parameter s, showing the direction of the source "electrical current" along 
the loop and the sense of the superflow encircling the loop. 
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O 

Figure 3. A two-dimensional vortex sub-manifold in a four-dimensional superfluid. Ti 
[= diR{s)] and T2 [— d2R{s)] are two coordinate tangent vectors to the sub-manifold 
at the point R{s). These tangent vectors define a "patch" of the sub-manifold. 




Figure 4. A segment of a vortex in a three-dimensional superfluid, showing the 
tangent vector dR/ds and two normal vectors and 7V^. 



Thus, for the case of a ^-dimensional vortex in a D-dimensional flow, the source 
takes the form 

Jd,-d^{x) = j d^Sea,...a^{daM ■ ■ ■ (da^Rd^) S^^'^X - R{S)), (11) 

where summations from 1 to are implied over repeated indices ai,a2, etc. We show 



m 



Appendix A that this form of the source does indeed describe a generalized vortex 



of unit vorticity, in the sense that the circulation of the superflow computed on any 
one-dimensional loop that surrounds the source a single time is unity (up to a sign that 
depends on the sense in which the loop is traversed). Furthermore, this form of the 



source obeys the condition of being closed, i.e., dd^Jd^. 



■dn- 



0, as we also show in 



Appendix A[ In the language of differential forms, Eq. ( JTTl) is expressed as 

j(x) = j d% 5^^\x - R{s)) (diR) A id2R) A ■ ■ ■ A (djpR). (12) 

It is also useful to consider the dual of the source, viz., the completely skew- 
symmetric rank-2 field 

^ Jdo-idoi-^) ■~ 'Jp^^^dxd2-djj,dD-\dD Jdx-dj^{x) 
^di--dr>-idr> 



j d%ea,...a^{daM " " " (da^Rd^) S^""^ {x - R{s)) . (13) 




Figure 5. A two-dimensional vortex sub-manifold in a four-dimensional superfluid. 
Shown are two tangent vectors, Ti and T2, and two normal vectors, and A^^, at 
the point R{s). 



Let US pause to analyze the factor 

from the integrand of Eq. (|T3l) . and hence to introduce some useful geometrical 
ideas. Being a skew-symmetrized (over the spatial indices di, d2 etc.) product of 
the set of vectors {daR}f=i tangent to the sub- manifold at the point R{s), the skew- 
symmetrization over indices associated with the sub-manifold coordinates (via ^ai - a^) 
is elementary, yielding a factor of so that the factor can equally well be written as 

^di--dD 

{d,R,,)---{djpR,^). (15) 

Being orthogonal to the plane tangent to the manifold (i.e., the plane spanned by the 
vectors {daR}f^i), this factor can be expressed in terms of a skew-symmetric product 
of any pair, and A^^, of unit vectors that are mutually orthogonal and orthogonal 
to the tangent plane (see Fig. H] for an example in D = 3 and Fig. [5] for an example in 
D = 4), i.e., as 

(Nl.Nl-NlNlJB, (16) 

in terms of a coefficient B. To determine B (up to a sign associated with orientation), 
we take the formula 

iNl_Nl - NlNlJB = e,,...,,idM ■ ■ ■ (dj^R,^), (17) 

and fully contract each side with itself to obtain 

B'^ = l^d^-dj^do-idD {diRdJ ■ ■ ■ (f^^-Rd^) ("^i-RjJ- ■ -{.djpR^^) 

= {^dxdi^d2d2 ■ ■ ■ ^d,pdjp ~ ^did2^d2di ■ ■ ■ ^djpdjp + ■ ■ ■) 

X {d,R,,y ■ ■{djpR,^){d,Rj^y ■ -{djpRs^) 
= + (OMidM {d2Rd2){d2Rd2) ■ ■ ■ idjpR,^)idjpR,^) 
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- {dMidM {dA,){d,R,,) ■ ■ ■ {djpR,^){djpR,J + ■■■ 
= detg, (18) 

where Qa'a are the components of the induced metric g on the vortex sub-manifold, also 
known as the first fundamental form [TT] , and are defined via 



dRdjs) dRdjs) 

Thus, for the geometric factor (fT4l) we obtain 

iNl_Nl - Nl^Nl), (20) 
and for the components of the dual source field -k J we obtain 

X {Nl.Nl - N^NlJ 6^^\x - R{s)) . (21) 



The measure d^ s A/det g is the infinitesimal area element on the vortex sub- manifold 
R{s). In the language of differential forms, Eq. f|T3|) is expressed as 

^ J(x) = j d%^d^{N^ A N^) 5(^)(a; - R{s)) . (22) 
3. Equilibrium superfiows of a given vorticial content 

In this section, we introduce energetic considerations and show that — provided we 
adopt a suitable choice of gauge for the rotational part V of the flow — the irrotational 
contribution to the flow V$ is zero in the equilibrium state. After that, introduce 
a convenient gauge potential A, and use it to solve for the flow in the presence of a 
given vorticial content, the latter being specified by a suitable source term. The gauge 
potential is a rank-0, skew-symmetric tensor, precisely analogous to the vector potential 
of three-dimensional magnetostatics. 



3.1. Minimization of superflow kinetic energy 

We take the energy of the system to be entirely associated with the kinetic energy of 
superflow. As we are addressing flows possessing a prescribed set of vortices, the energy 
associated with the absence of superfluidity in the cores of the vortices will not play a 
significant role. We assume that the mass-density of the superfluid p is constant, i.e., 
we take the fluid to be incompressible. The energy pE associated with the flow V can 
therefore be specified in terms of the functional 

E :=]- ! f/^xV,(x)Vrf(x), (23) 
^ Jo. 

which, for convenience, we shall refer to as the energy, despite the missing factor of p. 
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To determine the equilibrium flow in the presence of some prescribed vorticial 
content, we regard V in Eq. (jlj) as fixed and minimize E with respect to $. The 
corresponding stationarity condition reads: 

= tI^t = / d^x [V,{x) + V,$(x)] [V,{x) + VMx)] (24) 

= - Vd{VdHx)-Vd{x)), (25) 

where we have used the periodic boundary conditions on V and $ to eliminate boundary 
terms, and thus arrive at an overall condition of divergencelessness on the flow V. Now, 
V obeys Eq. ([6]), and is therefore not unique, being adjustable by any gradient of a 
scalar. This gauge freedom permits us to demand that the rotational part of the flow 
itself be divergenceless, i.e., 

VdVd = 0, (26) 

and we shall impose this demand. Thus, from Eq. 0251) we see that the minimizer $ 
is harmonic, i.e. V^Vd $ = 0, and subject to periodic boundary conditions. Thus, $ 
is a constant, and does not contribute to the equilibrium flow V, which is given by V 
alone, the latter obeying Eqs. IQ and fl26|) . Thus, the flow has energy (p/2) J d^xVdVd- 
(Minimization, and not just stationarity, of the energy follows directly from the positive 
semi-definiteness of the second variation oi E.) 

3.2. Fixing the complete flow: equilibrium state at given vorticial content 

To solve Eq. (I6l) in terms of the source and, more particularly, the location of the singular 
surface -R(s), we take advantage of the divergencelessness of V , Eq. fl26|) . This allows 
us to introduce a gauge field A, which is a completely skew-symmetric rank-^ tensor 
field, in terms of which V is given by 

Vdoi.^) = ^d^-dn ^dn-iM-djp{x). (27) 

We may subject this gauge field to a gauge condition, and the one that, for convenience, 
we adopt, is the transverse condition: 

Vd^Ad,...d^{x) = 0. (28) 

In the language of differential forms Eqs. ( l27|l and ( l28l) are expressed as 

V = ^dA and dA = 0. (29) 

To solve for A, we insert V, expressed in terms of A, into Eq. ([6]), to obtain 

(^di-do^dB-i(^d^-do-idD^dn-i^di-djp = ^Tl Jd^-dj^. (30) 

We then contract the two e symbols, which yields a [{D — 1)! termed] sum of products 
of [D — 1) Kronecker deltas, signed with the signature of the permutation tt, i.e., 

e^i-dD Qi-Jo-idD = X^- sgn(^) ■ ■ ■ (31) 
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and employ the skew symmetry of A and the gauge condition fl28p to arrive at the 
generahzed Poisson equation, viz., 

VdVdAd,...d^_^{x) = 27TJd,...d^{x). (32) 

By virtue of the adopted gauge condition, Eq. fl32l) has the cartesian components of A 
uncoupled from one another. Lastly, we introduce the D-dimensional Fourier transform 
and its inverse, defined as 



fiq) = / d^xe-'^-^fix), (33) 
Jn 

f{x)= [ d'^qe^'^-^fiq), (34) 



where g is a Z)-vector and d^q := d^q/ {2tt)^, and apply it to Eq. (!32|) to take advantage 
of the translational invariance of the Laplace operator. (We are implicitly taking the 
large- volume limit, and thus the inverse Fourier transform involves an integral over wave- 
vectors rather than a summation.) Thus, we arrive at an algebraic equation relating the 
Fourier transforms of the gauge field A and the source J, which we solve to obtain 

M-^A<i) = -^y-^^^^- (35) 

' ^^^^ Ip\ q-q ^ ' 

Inverting the Fourier transform, A{q), yields the real-space gauge field 

Ad,-dj^ix) = -j^ J d''q-—Jd,...d^{q), (36) 



from which, via Eq. fl27j) . we obtain the Fourier representation of the real-space velocity 
field, 

271- r „ iq^ e*^'^ 
Vdoi.^) = -^^d^-dn J A — -Jd^-dj^iq)- (37) 

For a source J given by a single vortex, as in Eq. (|TT1) . we have the following, more 
explicit form for the velocity field: 

2^ . f iD„ „-iq.x 



Vdo{x) = -^ ed,...dn J d ^^-ye 



j d^sea,...a,_,idaM ■ ■ ■ (da^Rd^) e^^-^(^^ (38) 



As shown in Appendix B the q integration may be performed, yielding the result 
^dD[x) 2 7r0/2 (^ - 1)! 

X J d^se,,..,^{daM ■ ■ ■ (da^Rd^) _ ^1^)1?" > (39) 

in which T{z) is the standard Gamma function 

Formula fl39l) for the velocity field V gives the structure of the superflow at any 
point X in the superfiuid in terms of the shape i? of a vortex. In the following sections, 
we shall deploy this formula to address the following issues: the velocity conferred by 
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this flow on vortices themselves, the velocity of hyper-spherical vortices, the dispersion 
relation for small deformations of hyper-planar vortices, the superflow kinetic energy 
for vortices of arbitrary shape, and the energy-momentum relation for hyper-spherical 
vortices. 



4. Dynamics of vortices in arbitrary dimensions 



4.I. Motion of vortices: asymptotics and geometry 

We have already established the flow field in the presence of a vortex, Eq. (!38l) . We 
now examine the velocity that this flow confers on points x = -R(S) on the vortex 
itself: 

1 r (0/2) 

x=i?(E) 27r^/2(0-l)!'''^-''^ 

(/?(S) - R{s)),^_ 



X 



■5 ^ai ■ 



2vr0/2 _ 

dajRdA-'^) ■ ■■9a, 



Rd 



(40) 



\R{T.)-R{s)f 

Observe that the integral in Eq. fHOj) is logarithmically divergent; physically, it is 
regularized via the short- distance cut-off ^. 

We now determine the leading contribution to the velocity in a regime in which 
the characteristic linear dimension of the vortex and the radii of curvature of the sub- 
manifold supporting the vortex are large, compared with ^. In this limit, the dominant 
contribution to the integral is associated with small values of the denominator. As we 



show in Appendix D the asymptotic behaviour of U{J1) is given by 

l_ T{Ip/2) 1 

m 

■ ■ (daj^Rdg,) idbdb>Rdr,_J gbb' 



UdJ^) 



X 



— {db'Rdo^i) (da^Rdi] 



^9bb' 



+ 2^dbRi) {dbd-yRd) {da.Rdr 



(dbda^RdJ ■■■{da^Rd^ 

■ (daj^Rdj^) {dyRdo-i) 

X {9bb9yy + 9by9bb' + 9bh'9w)] , (41) 

in which all terms, such as gradients of R, the metric tensor g and the inverse g of 
the metric tensor, are evaluated at the point S on the vortex sub-manifold, and the 
repeated indices b etc. are summed from 1 to Ip. The factor Ajp is the surface area of a 
^-dimensional hyper-sphere of unit radius (see, e.g., Ref. [12]). 

The asymptotic approximation to the velocity f/(S) of the vortex has a simple, 
geometrical structure to it. This structure becomes evident via the introduction, at each 
point S of the manifold, of a pair of vectors A^i(S) and A^^^S) that are normalized, 
orthogonal to one another, and orthogonal to the tangent plane to the vortex sub- 
manifold at -R(S), as made after Eq. f[T^ . The vectors A^^ and A^^ span the manifold 
of dimension two complementary to vortex sub-manifold, playing the role played by the 
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normal vector of a hyper-surface (i.e., a sub- manifold of codimension 1). The choice of 
the vectors A^^(S) and A^^(S) is not unique; alternative linear combinations of them, 
provided the pair remains orthonormal to one another and orthogonal to the tangent 
plane, will do equally well. However, what is independent of this choice, is unique and, 
from a geometric perspective, is the natural object, is the skew-symmetric product of 
iVi(S) and iV2(E), viz., the two-form 

N := A N^, (42) 

the components of which are given by 

N,,,, := Nl Nl - Nl Nl , (43) 



already featuring in Eq. f|T6l) et seq. As we show in Appendix E, the right-hand side of 
Eq. ( HTl) can be expressed in terms of the components of the A^^ A A^^ two-form, which 
then reads 

C/d(E) ^ ^ Qaa' N,,> {dada'Rd') HL / ■ (44) 



Note that we have restored to dimensional form of the velocity in Eq. fl44j) . according 
to the prescription given in the footnote 

Observe the role played by the geometry here. The local shape of the vortex enters 
through the intrinsic geometric characteristic, the metric g [i.e., the first fundamental 
form, Eq. (IT9]) ]. as well as through the extrinsic geometric characteristic Nm' dada'Rd' 
(i.e., the second fundamental form [H]). For the familiar case of D = 3, one choice 
for the two vectors A^^ and A^^ is the normal and the binormal to the vortex curve. 
Together with the tangent to the 1-dimensional vortex line they form an orthonormal 
basis (i.e., a Frenet-Serret frame), which reflects the extrinsic geometry of the vortex 
line. In contrast with its higher- dimensional extensions, the vortex line does not possess 
any intrinsic geometry. 

4-2. Velocity of circular, spherical and hyper- spherical vortices: exact and asymptotic 
results 

We now examine the velocity of vortices for the illustrative cases of circular, spherical 
and hyper-spherical ^-dimensional vortices of radius L, embedded in Z)-dimensional 
space. In such cases, the vortex sub-manifolds can be taken to have the following form: 

R{s) = L(sin si sin S2 ■ ■ ■ sin S0_i sin s^, 
sin si sin S2 ■ ■ ■ sin Sjp^i cos Sjp, 
sin si sin S2 ■ ■ ■ cos Sjp^i, 



sm si cos S2, 
cos Si, 



§ To restore the physical units, we multiply velocities by h/M, energies by pfi^ /A'P, momenta by ph/M 
and frequencies by h/M. 
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0), (45) 

with the Ip parameters (si,...,^^) respectively ranging over intervals of length 
(tt, 27r, . . . , 27r), and where, for each D, we have chosen the manifold to be centered 
at the coordinate origin and oriented so that it is confined to the xn = hyper-plane. 
For example, in the case of D = 3, the vortex sub-manifold is then a ring of radius 
L, lying in the X1-X2 plane. As the vortices are circular, spherical or hyper-spherical, 
rotational invariance ensures that all points on them move with a common velocity, and 
that this velocity points along the direction perpendicular to the [D — l)-dimensional 
plane in which the vortices are confined, i.e., the D direction. 

Formally, we can write an exact expression for the velocity f/(S) of any point S 
on the manifold, by substituting R{s) into Eq. pUj) . Thus, may choose any point on 
the manifold at which to compute the velocity, and for convenience we choose the point 
given by si = [i.e., R = (0, . . . , 0, L, 0)]. Two convenient vectors orthogonal to each 
other and normal to the vortex sub-manifold there are: 

A^^ = (0,0,..., 0,0,1) and A^^ = (0, 0, . . . , 0, 1, 0). (46) 

By substituting R, and A^^ into Eq. (1411) we obtain 

Uni^) = ^ln{L/0, (47) 
regardless of the dimension D. 



4-3. Motion of weak distortions of linear, planar and hyper-planar vortices 

Next, we assume that the manifold is a weak distortion of a ^-dimensional hyperplane, 
and examine the motion of this distortion that results from it being transported by the 
flow. Such hj^er-planar vortices are the generalization to the D-dimensional setting of 
the one-dimensional straight-line vortex in the three-dimensional setting — flrst studied 
for normal fluids by Thomson in 1880 |T5] — and the flat plane vortex in the four- 
dimensional setting. The points on the undistorted hyper-plane are given by 

X = F{s), (48) 

where F{s) = J2f=i Sa^a, the set {ed}d=i is a complete orthonormal Euclidian basis in 
the ambient space, and the variables s range without bound. To obtain a weak distortion 
of the flat manifold we augment F{s) with a "height" function, 

H{s) = Hd-i{s) eo-i + Hd{s) eo, (49) 

which has nonzero components i/z)-i and Hu only, i.e., it points in the two directions 
transverse to the flat manifold. Figure [6] illustrates the weak distortion of a plane in 
D = A. Thus, the points on the weak distortion of the flat manifold are given by 

X = F{s) + H{s). (50) 

We now analyze the dynamics of the small perturbations (i.e., ripples) that 
propagate on the manifold, by constructing an equation of motion for them. The flrst 
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Figure 6. A weak distortion of a hyper-planar vortex in a four-dimensional superfluid. 
The height function H characterizes the distortion, which is perpendicular to the 
undistorted plane. 



ingredient is the velocity U of any point a on the manifold, which we determine by 
applying Eq. ( l40l) . subject to the replacement 

R{s) ^ F{s) + H{s). (51) 

Thus, we arrive at the flow velocity at the point x = R{a): 



r(i/)/2) 

t/,,(a) = 27r^-5(^e,,... 



sen 



X da, [F,, is) + Ha, (s)) ■ ■ ■ da^ {F,^ is) + H,^ 
{F{a)-F{s) + H{a)-H{s)),^_^ 



») 



\F{a)-F{s) + H{a)-H{sr ^^^^ 
We are concerned with small distortions, and we thus expand the right hand side 
of Eq. fl52l) to first order in the height H; the details of this expansion are given 
in Appendix F For all points cr of the manifold, the resulting velocity U{a) points 



in a direction transverse to the flat manifold, its only nonzero components being 
approximately given by 



«.r(^/2) 



\F(a) 



X 



-F(s)r 

[H{a)-H{s))^-J2 



a 



OH^js) 
dsu 



(53) 



where the indices 7 and 7' (which, are summed over when repeated) take only the values 
[D — 1) and and e-^y is the skew-symmetric tensor with e^^i^o := -|-1 

We pause to note that the superflow velocity field V{x) is translationally covariant, 
in the sense that a global shift in the position of the vortex sub-manifold by a 
constant vector C (i.e., R ^ R + C) produces the shift of the superflow velocity field 
V{x) V{x — C). The vortex velocity given by Eq. (l52l) . and its approximation 
Eq. (1531) . both transform in accordance with this requirement. In particular, it should 
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be noted that Eq. fl53p is invariant under translations of the distortion of the form 
H^{a) = H^{a) + C^. 

Returning to the calculation of the vortex velocity for the case of weak distortions 
from hyper-planarity, the form of the denominator in the integral in Eq. (!53l) suggests 
that the dominant contribution to f/-y'(S) comes from short distances, i.e., s ~ a. We 
therefore approximate this integral by making a Taylor expansion of the term H{a) 
about cr = s to second order, thus arriving at the approximate vortex velocity 

,r(^/2) f d^s 



x(.-.Ua-.)4^. (54) 

To obtain a closed system of equations for the dynamics of the height function, we 
observe that the superflow carries the vortex sub-manifold with it, so that the evolution 
in time t of the height is given by 

dH^{a,t)/dt = U^{a,t). (55) 

In the approximation (1541) . this equation of motion becomes a first-order differential (in 
time), nonlocal (in "internal space" i.e., in the coordinates a that span the undistorted, 
hyper-planar vortex sub-manifold) system of two coupled linear equations that is 
translationally invariant in both (internal) space and time. As such, it is reducible 
to a two-by-two matrix algebraic form under Fourier transformation in both (internal) 
space and time, which we define as follows: 

L{q,uj) = [ d^a e"^'?-" / dt e''^' L{a,t), (56) 



Computing the spatial Fourier transform of Eq. ( 1551) is addressed in Appendix G the 



temporal Fourier transform is elementary. Introducing, for convenience, the coefficient 



M\k\)--=Y^^k'\n{l/m, (58) 



arrived at via the spatial Fourier transform (see Appendix G ), the equation of motion 
becomes 

— iujH^{q,uj) = A(|g|) e^y iJy (g, cj). (59) 
The solvability condition for this system is 

det ( Z n ^^'''^ ] = 0' (60) 

and gives rise to the following dispersion relation for the spectrum of normal modes of 
oscillation, with two modes for each value of the ^-vector q conjugate to the position 
in the fiat vortex sub-manifold hyper-plane: 
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where we have restored the dimensional factors. The normal modes are themselves given 
by 



^^[<l-T^t)^ (62) 



and these correspond to left and right circularly polarized oscillations. 

5. Superflow kinetic energy for given vorticial content 
5.1. Exact evaluation 

We now use the solution for the velocity field, Eq. (1371) . together with the equation for 
the kinetic energy E of the velocity field, Eq. (123|1 . to obtain the kinetic energy in terms 
of the source J: 

^ ^ V,{x)V,{x) = \ 



E=^ I d^'x VJx) VJx) = ^ d^'x V:(x) V, (x) 



q' ■ q' 



jfxe,,...,^jd-q'q',^_^e-^'- 
J Q-Q 



2 



X / d^q J},..d'^iq) Jd.-d,{q) ^^^t|?^, (63) 

where 5^^\q) := {2-1^)^ 5^^\q) . For a source J given by a single vortex, as in Eq. ( fTTl) . 
we then have the more explicit form for the energy: 

X j d^Sea,...a^{da,RdAs))---{da^Rdjj,{s)) 

Writing X for (-R(s') — R{s)) and using the result for the q integration given in 



Appendix C, the energy becomes 

2 



^<i[-d'^ed,...djd'^dn J d^s'£a[-a'^ida'^Rd[is')) ■ ■ ■ (9„^i?rf^ (s')) 
X j d^S ea,...a^ida,RdAs)) ■ ■ ■ (da^Rd^is)) 
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{X ■ x f" 

where Vi and V2 depend only on D and are given in Appendix C 



{S.^^.^o-. V, + \\ P.) , (65) 



The final step in constructing the energy of a single vortex involves contracting the 
D-dimensional e symbols and taking advantage of the skew-symmetry of the sources. 
Thus, we arrive at the result 

X I d^s'^.;.,,,., (4; fi*(s'))---(4i,_.fi.,_,(.5'))(9«i,-Rji, (*>')) 
X J <i*s£„,.„„,(3„,fij,(s))---(3„,_,fij,_,(s))(9„,fij,(s)) 

The factor Ad-i is the surface area oi a {D — l)-dimensional hyper-sphere of unit radius 
(see, e.g., Ref. [12]), and T{z) is the standard Gamma function (see, e.g., Ref. [H]). 

Formula ( |66ll gives the kinetic energy in terms of the local geometry of "patches" 
of the vortex sub-manifold and interactions between such "patches'.' The asymptotically 
dominant contributions to the energy come from nearby patches, and we evaluate these 
contributions in the following section. The linearity of the theory ensures that the 
energy associated with configurations involving multiple, rather than just one, vortices 
is obtained via an elementary extension, comprising self-interactions terms, such as the 
one just given, for each vortex, as well as mutual interaction terms for pairs of distinct 
vortices. 



5.2. Asymptotics and geometry 



We now analyze the energy of a vortex asymptotically, subject to the conditions 
discussed earlier in Section 14.11 (i.e., the characteristic linear dimensions and radii of 
curvature of the vortex are much larger than the short-distance cut-off C,)- In this 
regime, the dominant contribution to the energy is associated with the large values 
of the integrand that occur when the (^-dimensional) integration variables s and s' 
are nearly coincident, so that the denominator (X ■ X)^/^ is small. To collect such 
contributions, it is useful to change variables from {s, s') to (S, a) as follows: 

Sa ■= Sa, (67) 

<:=S„ + a„ (68) 
the Jacobian for this transformation being unity, 
d{s,s') 



det 



a(S,a) 

so that Jd^s' Jd^s... = Jd^Jd^a... . 



(69) 
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The energy formula fl66l) comprises two pieces, and we now consider them in turn. 
The first, (a), involves the integral 

/rf^s I ^^"^^^^ e,^...,^ {d,R,,{j: -.))••• {d,^R.,{j: - s)) 

Xea,...a^{daMs))---{da^Rd^{s)). (70) 

The second, (b), involves the integral 



X ea,-a^{da,RdAs))---{da^_,Rd^_,{s)){da^Rd^{s)). (71) 

Our strategy will be to Taylor-expand the integrands in terms (a) and (b) about the 
point (7 = 0, keeping only the leading-order terms. 

Let us focus on term (a). It is adequate to approximate the numerator at zeroth 
order: 

e,....,.^(9,. /?,,(S))- • id,.^Rd^{T.)) e,,...,^ {d^M^))- " ■(9,^i?d^(S)) . 

Note the occurrence of Ip factors of the induced metric g^^i, defined in Eq. (1191) . In 
terms of the numerator of term (a) becomes 

and, by the well-known formula for determinants, this is 

det g . (73) 

To complete the approximate computation of term (a) we examine the denominator, 

(X-X)"^/' = |i?(S + (T) -i?(S)r^. (74) 

By Taylor-expanding for small a and retaining only the leading-order contribution, we 
find that the denominator becomes 

{(^a' ga'a(^ay^^^ ■ (75) 

Thus, for term (a) we have the approximate result 

. (76) 

A convenient way to analyze the integration over a is to make a ^-dimensional 
orthogonal transformation from a to coordinates C that diagonalize g, the eigenvalues 
of which we denote by {7a}f=i, and then to "squash" these coordinates to C, (defined via 
Ca := JlaCa)) SO as to "isotropizc" the quadratic form in the denominator. Following 
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these steps, postponing to the following paragraph a discussion of the integration limits, 
and omitting an overall factor of our approximation for term (a) becomes 

rf^Edet5(S) ^ '^^'^ 



(cTa' 9a' a{^) CTa)^^^ 

^ J ci^Sdet5(S) j- "^^^ 



Ef=l7a(S)a 

det5(S) f d^C 



0/2 



= j d^E Vdet^(E)^0 j 



Z_/a=l So 

dzz"^-^ 



0/2 



where Z := \jY.a=iQ- 

Up to now, we have been vague about the limits on the integrations over the 
(7 variables. In the short-distance regime there is, for each value of E, a physically 
motivated cutoff that eliminates contributions from regions of a in which the separation 
|i?(E + (t) — -R(E)| is of order ^ or smaller. This cutoff is associated with the fact that 
our overall description is not intended (and is indeed unable) to capture physics on 
length-scales associated with the linear dimension ^ of the core of the vortices. In the 
long-distance regime, there are integration limits associated with the finite size of the 
vortex, which we take to be of order L. (Thus, we rule out of consideration vortices 
whose shape is strongly anisotropic; already ruled out are, e.g., rough vortices, for 
which the Taylor expansion of the functions describing the shapes of the vortices would 
be unjustified.) In particular, the last step in our approximate computation of term (a), 
in which we transformed to the single radial coordinate in the ^-dimensional fiat plane 
tangent to the vortex at the point s = E, is only valid if the limits of integration possess 
the rotational symmetry of this tangent plane, which they typically do not. (The case 
of D = 3 is an exception.) However, to within the logarithmic accuracy to which we 
are working, it is adequate to ignore these complications and approximate what has 
become the Z integral by ln(L/,^). Any refinement would just change the argument of 
the logarithm by a multiplicative numerical factor, and this would produce an additive 
correction to our result that is of sub-leading order. Thus, our approximation to term (a) 
becomes 

U)\ j d^EVdetg(E)^0ln(L/e). (77) 

The remaining integral, / d^E •\/detg(E), is famifiar from differential geometry, and 
simply gives the ^-dimensional volume of the sub-manifold supporting the vortex, which 
we shall denote by [and which is a functional of R{-)]. Thus, our approximation to 
term (a) becomes 

Ip\AjpWjpHLli)- (78) 
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A similar procedure allows us to find an adequate approximation to term (b), 
as we now show. In contrast with term (a), owing to the two factors of '■ = 
+ cr) — the numerator vanishes at zeroth order in a, and requires Taylor 

expansion, 

Xd^aadRd/dJla, (79) 

to identify the leading-order contribution to the numerator. The remaining factors in 
the numerator can be evaluated at zeroth order, which then becomes 

^a'l-a^ da[Rd^ ■ " " da'^_^Rd^_, ^a'^Rd'^^ ^a^-a^ da.Rd, ' " " da^_,Rd^_, ^a^j^Rd^ 



= (0 - 1)! Saj^a'^_^ Cr^'^_^ 9a^aD-i ^a^-i = ^- ^^a' 9a'a (^a , (80) 

in which all terms involving R and g are evaluated at the point S, and where we have 
made use of the formula for the determinant, 

^a'^-a^ Qa'^ai S'a^aj " " " da'^^ajp = ^^ai - a^ det g, (81) 

and, once again, introduced the induced metric ga'ai^) as well as a contraction 
^ai---a0_j^a£ai---a0_-^a' = (-0 — ^V-^aa'- As for the the denominator of term (b), this we 
Taylor-expand in the manner already used for term (a), noting that the overall powers 
in these denominators differ by unity. Thus, apart from an overall factor of (-0 — 1)!, 
we have for term (b) the approximate result 

^a' 9 a' a 



rf^Sdet^(S) j d^a 



d^Sdet^7(S) / (82) 



{pa' 5'a'a(S)(Ta)^/^ 

which is of precisely the structure that arose in the analysis of term (a), and thus the 
approximate result for term (b) becomes 

(0-l)!^0W0ln(L/O. (83) 

Finally, we assemble terms (a) and (b) to arrive at an asymptotic formula for the 
energy of the equilibrium flow in the presence of a vorticial defect, which reads 

E ^ ^v^r (f - i) r (f - 1) An-x Ajp ln(L/0 W0. (84) 
This formula simplifies, using the equation for Ref. 112]), to give 

(p]^)^>V0ln(L/O, (85) 
where we have restored the dimensional factors. 
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5. 3. Elementary argument for the energy 

In the previous subsection, 15.21 we have shown that, in the asymptotic hmit of a large, 
smooth, ^-dimensional vortex sub-manifold having the topology of the surface of a 
{D — l)-dimensional hyper-sphere, the kinetic energy depends on the geometry of the 
sub- manifold only through its ^-dimensional volume Wjp. This result is consistent with 
the following elementary argument. 

Imagine starting with such a vortex sub-manifold, and "morphing" its shape into 
that of a hyper-sphere having the same ^-dimensional volume. Next, cut the hyper- 
sphere along an equator, so that it becomes two hyper-hemispheres of equal size. Then, 
distort the hyper- hemispheres into flat hyper-planes separated by aL (i.e., a distance 
of order L), maintaining their ^-dimensional volumes at ^Wjp, and imposing periodic 
boundary conditions on opposing sides of the hyper-planes. This last step is bound 
to involve some stretching of the sub-manifolds, but let us assume that this will only 
affect the energy sub- dominant ly. We choose the separation of the hyper-planes to be 
0{L), so as to best maintain the geometry of the original sub-manifold. We can think 
of the resulting D-dimensional system as comprising a ^-dimensional stack of identical, 
two-dimensional, films of superfiuid having thicknesses ^ in every one of their Ip thin 
dimensions, each film containing a point-like vortex-anti-vortex pair. 

Now, owing to the translational and inversion symmetry of the stack, the flow in 
each film is confined to the film and is simply that associated with vortices in the film, 
and therefore the energy of the flow is given by the formula: 2tt^^ ln(aL/^) (see, e.g., 
Ref. p]). Furthermore, the number of films is ^W^/^^. Thus, consistent with Eq. ( l85l) . 
the total energy is given by 

^x2,f*l„(i^)=.W^|l„| + 0(l)}, (86) 

where the 0{1) correction is associated with the factor a. 

Sub-leading corrections to this result for the energy will not depend solely on the 
vortex sub-manifold volume but also on its more refined local and global geometrical 
characteristics. 



5.4- Energy of circular, spherical and hyper- spherical vortices: exact and asymptotic 
results 

In this section, we analyze the energy of a geometrically (and not just topologically) 
hyper-spherical vortex sub-manifold. The function R{s) for such a vortex is given by 
Eq. psl) . and by substituting this form into Eq. (165!) we obtain a formally exact result 
for the energy. 

To illustrate this point, consider the case of a two-dimensional vortex in D = 4, for 
which the energy reads 

E = -L? I dsi ds[ / ds2 ds^ cos(s2) cos(s2) 

O J-K J—k/2 
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/ 2 

\ 1 — cos(si — s[) COS S2 COS ^2 — sin S2 sin 

As with the case of arbitrary D, the integrand in Eq. flHTl) is singular for (si, S2) = {s'l, S2); 
therefore, we follow the procedure, described in detail in Section 15.21 of expanding the 
integrand around this point of singularity, hence obtaining the asymptotic result that 
the energy is given by 

E ^ An^ L"^ \n{L/^). (88) 

Turning now to arbitrary dimension D, the energy of a hyper-spherical vortex can 
be obtained by replacing the volume of the manifold in Eq. fl85l) by the surface area of 
a {D — l)-dimensional hyper-sphere, viz., Ad~iL^ (see Ref. [12]), so that we arrive at 

E^7iAd-iL^ HL/^. (89) 




5.5. Energy-momentum relation for vortices 

We now use the results for the energies and velocities of circular, spherical and hyper- 
spherical vortices, obtained respectively in Secs. l4.2 land l5.4[ to construct the momentum 
P of such vortices and, hence, their energy-momentum relations. To do this, we observe 
that the energy E{L) and velocity V{L) depend parametrically on the radius of the 
vortex L. Hence, we may construct the parametric dependence of the momentum P{L) 
via the Hamiltonian equation 



V{L) = ^E{p) 
in view of which we have 



P=P{L) 



(90) 

dP/dL' ^ ' 



where we have retained only the leading behavior at large L. By integrating with respect 
to L we obtain the parametric dependence of the momentum on the vortex radius, 

PiL)^2n^^An-iL''-\ (92) 

and by eliminating L in favour of E we arrive at the energy-momentum relation 

E^ri{D)P^ [In^ + ln^] , (93) 



where r]{D) is a D-dependent constant given by 

V{D)=^^^^^^^ (-^ , (94) 



D-l \2D-A^ 

and the sub-leading logarithmic term is sensitive to the number of particles per coherence 
volume. 



Quantized vortices and superflow in arbitrary dimensions 



25 



6. Scaling and dimensional analysis 

To ease the presentation, we have chosen to drop factors of h, M and p from our 
derivations, restoring physical units only when giving results at the ends of the various 
sections. It is worth noting, however, that many of our results can be obtained via 
scaling and dimensional analysis, up to overall dimension-dependent factors and cutoff- 
dependent logarithms, as we now discuss. 

There are two characteristic length-scales: the linear dimension of a vortex L, 
and the vortex core size ^. The latter serves as a cutoff, and therefore only features 
in logarithmic factors, which are inaccessible to our scaling and dimensional analysis. 
There are two characteristic mass-scales: the mass of each of the condensing particles M, 
and pL^, where p is the mass-density of the superfiuid. Lastly, there is a characteristic 
frequency-scale: h/ML'^. 

To ascertain the scaling form of the velocity of a vortex, we note that its derivation 
involves only the generalized D-dimensional Ampere- Maxwell law, and thus is insensitive 
to p. Hence, dimensional analysis suggests that 

V Lx — ^ —, (95) 

ML2 ML' ^ ^ 

which agrees, e.g., with Eq. (jlTj), except for numerical and logarithmic factors. As for 

the energy of a flow associated with a vortex, it is proportional to the density p, and 

therefore dimensional analysis suggests that 

which agrees, e.g., with Eq. fl85l) . except for numerical and logarithmic factors. 
Regarding the momentum, it was computed from the Hamiltonian equation of 
motion (pOl) . and thus, like the energy, scales as p^. Dimensional analysis therefore 
suggests that 

P.,i«xLx-A_~,L-A, (97) 

which agrees, e.g., with Eq. ( l92i) . except for numerical and logarithmic factors. 
Concerning the frequency of oscillations of distortions of wave-vector g of a hyper- 
planar vortex [see Eq. (EH)], these are driven by the flow velocity, and therefore do not 
depend on p. Moreover, the undistorted vortex does not have a system-size-independent 
characteristic size; instead, the characteristic length-scale is set by q~^. These scaling 
notions, together with dimensional analysis, suggest that 

i^-^q', (98) 
consistent with Eq. (|6T|) . except for numerical and logarithmic factors. 
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7. Concluding remarks 

We have explored the extension, to arbitrary spatial dimension, of Onsager-Feynman 
quantized vortices in the flow of superfluid helium-four, focusing on the structure and 
energetics of the associated superflow and the corresponding dynamics of the vortices. 
To do this, we have analyzed the superflow that surrounds the vortices at equilibrium 
by invoking an extension, to arbitrary dimensions, of the three-dimensional analogy 
between: (i) the magnetic fleld of Ampere-Maxwell magnetostatics in the presence of 
specified, conserved electric currents concentrated on infinitesimally thin loops; and 
(ii) the velocity field of equilibrium superfiuidity of helium-four in the presence of 
specified, conserved vorticity concentrated on infinitesimally thin vortices. 

By constructing the appropriate conditions that the velocity field must obey at 
equilibrium, if the flow is to have a given vorticial content, we have determined the 
corresponding flow in terms of this vorticial content, via the introduction of a suitable, 
skew-symmetric gauge potential of rank two fewer than the spatial dimension. The 
structure of the flows associated with generically shaped vortices has enabled us to 
develop results for the dynamics that vortices inherit from these flows in all dimensions, 
and for the way in which this dynamics rcficcts both the intrinsic and extrinsic geometry 
of the vortices when the vortices arc large and smooth, relative to the vortex core size. 
This structure has also allowed us to ascertain the velocities with which the higher- 
dimensional extensions of circular vortex rings (i.e., vortex spheres and hyper-spheres) 
propagate, as well as to analyze the higher-dimensional extension (to planes and hyper- 
planes) of the Kelvin problem of the vibrations of a line vortex. 

For all dimensions, we have expressed the equilibrium kinetic energy of superflow 
in the presence of a generic vortex in terms of the shape of the vortex, and have 
identified the leading contribution to this energy — for vortices that are smooth and 
large, relative to the vortex core size — as being proportional to the volume of the sub- 
manifold supporting the vortex. We have also computed the kinetic energy for the 
higher-dimensional extensions of vortex rings, and used these energies to determine the 
momenta, and energy-momentum relations, of these maximally symmetrical vortices. 
In addition, we have shown how the structure of our results, if not their details, can be 
obtained via elementary arguments, including scaling analyses. 

The geometry of the vortices — both the intrinsic geometry of their shapes and 
the extrinsic geometry of the manner in which the vortices are embedded in the high- 
dimension ambient space — plays a central role in our developments. The natural 
language for such developments is that of exterior calculus and differential forms, but 
we have not made this language essential, using instead the equivalent, but more widely 
known, language of skew-symmetric tensor fields. 

Further issues that one might consider within the present framework include lattices 
of higher-dimensional extensions of vortices associated with equilibrium states of a 
rotating superfluid, along with the collective vibrations of such lattices (i.e., Tkachenko 
modes), as well as Magnus-type forces on vortices provided by background flows. One 
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might also consider the energetics and motion of knotted vortices in three dimensions 
(cf. Kelvin's "vortex atoms" [16]) and — where mathematically available — topologically 
more exotic, higher- dimensional vortices, including higher-genus surfaces. Moreover, 
one might apply the the energetics of higher-dimensional extensions of vortices, and 
the interactions between them, to a re-examination of the statistical mechanics of 
the superfluid-to-normal phase transition (i.e., the "lambda" transition) in higher 
dimensions. Still more challenging would be the generalization of the circle of ideas 
discussed in this paper to the non-linear setting provided by non-Abehan gauge fields. 
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Appendix A. Source terms and quantized circulation 

In this appendix, we discuss two properties of the source field J. First, we show that the 
circulation, Eq. ([21), of the superflow around the vortex sub-manifold associated with 
the source term given in Eq. ffTTj) is quantized to unit value. We then show that the 
source field is divergenceless. 




Figure Al. A closed path dAi around a segment of a vortex loop specified by R{s), 
with arc-length parameter s. dM. is the boundary of the two-dimensional region Al, 
which is punctured by the _0-dimensional vortex. 

We begin by considering a closed, one-dimensional path dM., centered on the vortex 
sub- manifold and encircling it, along with the associated two-dimensional region M. 
bounded by DM. and is punctured by the vortex sub-manifold, as shown in Fig. lAll In 
the language of differential forms. Stokes' theorem, applied to Eq. ([9]), determines that 
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the circulation k around dAi is given by 



K= (b Vd{x)dxd= <h V (A.l) 

JdM JdM 

dV = 2tx I icj=^ {-k J)did2{x) dxd^ A dxd2 ■ (A. 2) 



M JM 2! 

We take dAi to be circular, with the points on it being given by 

X = R{so) + aiN^ + a2N^, (A.3) 

where ai and 02 vary subject to the constraint + a| = const. > 0, and A^^ and 
A^^ are any pair of mutually perpendicular unit vectors that are perpendicular to the 
vortex sub- manifold at the point Sq. An elementary computation then shows that 
dxdi A dxd2 = dai A da2 {N^ A N'^)d^d2i so that, from Eq. (lA.2p . the circulation becomes 
27r f 

ft; = ITT / (* J)d,d2 {x) (N^ A N^)d^d2 dai A rfag. (A.4) 
^- Jm 

Next, we insert the explicit components '^Jdid2i Eq. f|T5]) . for the specific form of M. 
defined by Eq. (lA.Sp . to arrive at 
27r 



K - 

2 



\ j^j d'^Sed^.-.d^dn-^dAdaM ■ ■ ■ {da^Rd^) {N^ A N^)dr,^,dr, 

X 5^^) (i?(so) + aiN^ + aaA^^ - R{s)) dai A da2. (A.5) 

Next, we Taylor-expand R{s) in the argument of the delta function about the point 
s = So and use the orthogonality of A^^ and A^^ with the tangent vectors {dR/dsa}f=i 
at So to obtain 

2n 12 

K = —edi-dj^dD^ido{9aiRdJ ■ ■ ■ (daj^Rdjf,) {NAN )do-ido 



j d'^s5^'>-^^{{s-So)adaR{s)\s=. 

X ! 5^^\ai)dai ! 5^^\a2)da2. (A.6) 



2! 

X 



The last two integrals each give a factor of unity. The remaining one, over s, is readily 
seen to give a factor of 1/| det(9ai?)| provided we invoke the fact that the term daR is 
a square Ij) ^ Jj) matrix in the subspace of vectors tangent to the vortex sub-manifold. 
As we discuss in Section [221 the factor ed^-d^dD-idoidaiRdi) ■ ■ ■ (daj^Rd,^) is given by 
(A^^ A N^)dD-idD det{daR). Hence, in the formula for the circulation the determinants 
cancel, and the remaining two factors of the form A^^ A A^^ contract with one another 
to give a factor of 2, so that we arrive at the result 

K = ±1, (A.7) 

the sign depending on the sense in which the encircling path dAi is followed. Thus, 
we see that the circulation of the superfiow around the vortex sub-manifold associated 
with the source term given in Eq. ( ITTl) is indeed quantized to unit value. 

In order to discuss the divergencelessness property of J, we should distinguish 
between two cases. In the case of D = 3 the source field is a vector, which consequently 



Quantized vortices and superflow in arbitrary dimensions 29 

has no skew-symmetric properties. In arbitrary dimension D greater than or equal to 
4, the source field J has a tensorial structure that we exploit in our analysis. Let us 
consider first the case case oi D = 3. The source term is then given by 

Mx) = I ds^^ I ^e--(-«W), (A.8) 

and its divergence may be manipulated as follows: 



isj 


d^q . 


(2vr)3^ "^'^ 


d^q 


[ ds (e 


(27r)3 


J as 



ds 



(A.9) 



The integrand in the last line of Eq. (1A.9I) is a total derivative of a periodic function of 



s, and therefore the integral vanishes. 

For D > A the source term has the structure 

1D 



and its divergence may be manipulated as follows: 



j d^SEa^-ajfida^Rdi) ■ ■ ■ (daj^Rdjp, 



d^q 



(27r)^ 

- e-«-(--«(^))e.,...,^9,J(9,,i?,J ■ ■ ■ {da.Rd,)]}- (A.ll) 
In the last line of Eq. (lA.lip the first term vanishes because the integrand is a total 



derivative of a periodic function, whilst the second term vanishes via the skew-symmetric 
properties of Levi-Civita symbol. Hence, we have the divergencelessness condition 
^diJdi-dg, = 0. 

Appendix B. Computing the velocity field 

As shown in Section [221 the velocity field due to a single, ^-dimensional vortex is given 
by 

27r f iqdn-^ ig.^ 



'sea,...a,_,id.M ■ ■ ■ (daM e^^'^^^^. (B.l) 
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Performing the q integration using Eq. (lC.2p of Appendix C, and noting the formula for 
the surface area of a D-dimensional sphere of unit radius ^2] , the velocity field becomes 



2tt 



J Q-Q 



X ddi,_, I d^q 



Jq-(R(s)-x) 

q-q 



712 



-D 



-Ad-1 r(f - 1) r(-|) ed^...do 



X j d^Sea,-a^ida,RdJ ■ ■ ■ (da^Rd^) (Srf^. Jx - i?(s) | ^) 

^ (x — R{s))^^__^ 

1 r(^/2) 

X J d^Sea,...a^{daM ■ ■ ■ (da^Rd^) _ • (B-2) 

Appendix C. Invariant integrals 

Consider the q integral in Eq. (!64|) . Writing X for {R{s') — R{s)), and invoking 
rotational symmetry and dimensional analysis, we have 

qd 1 



(27r) 



Od'd Vl + V2 



Contracting first with 6 and then Xd> Xd, we arrive at the invariant integrals 



D^JqX 



d'^qe 



1 



g-g ~ {X-Xf^^'' 



(C.l) 



(C.2) 



(C.3) 



[X-X)]-^^^ {q.qf (X-Xf/'' 
The integrals flC.2p and flC.2p can be computed in hyper-spherical coordinates, and we 
thus obtain for Vi and V2'- 

v,= ^z.-iv^r(f -i) 2^-^r(f -1), (C.4) 

= (2 - D) Ad-1 v^r (f - i) 2"-' T (f - l) , (C.5) 

where Ad is the surface area of a D-dimensional sphere of unit radius [12], and T{z) is 
the standard Gamma function 
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To determine the asymptotic behaviour of the velocity, we start with Eq. (I40l) . which 
contains the factor 



Sea,-aada,RdAs)---daaRdJs) 



and Taylor-expand the integrand about the point S to obtain 

(daiRdi) ■ ■ ■ {dajpRdjp)\s=j: 
+ {daM ■ ■ ■ {dbda^RdJ ■ ■ ■ {da^Rd^)l^A^t " Sfe) + 
-(<9fe^do-i)L=E(^b-^b) 



- -{dbdb'Rd^_^)\^^j.{sb - - St,) + 

- idbR^)U^{sb - J^b) 

- ^{dbdb'Rd)\s=A^b - ^b){sb' - Sft/) H 

- {dbRd)l=Asb - ^b) 

- l:{dbdb'Rd)\,=j^{sb - T.b){sb' - Sfc/) H 



(D.l) 



X 



X 



D 

" 2 



(D.2) 



Next, we expand the terms in square brackets, retaining only terms to leading order, 
and collect terms of similar structure in the integration variables s, observing that two 
such structures arise: 



01020314 



Sai)(s 



12 ^0-2 ) 



i^ai Sc[-^)(Sq,2 ^a2){^a3 ^as)^^^ ^04; 



The asymptotic behaviours of these structures are given by 



(D.3) 
(D.4) 

(D.5) 



7"(4) (\^\ ~ "^0 9aia2 Qajai ~^ Qaias 9a2a-A ~l~ 9aia4 9a2a3 ]r.(T I C\ (T^ P.\ 

where g := g~^, and g and g are evaluated at S. 

Appendix E. Geometric content of the vortex velocity 

Equation fl^Tl) for the asymptotic approximation to the vortex velocity ^/^^(S) contains 
three terms, 



•^dn '■- ^d^-do^ai-aj.(,da^Rdi) " " " (Oa^Rdj^) {dbdb' Rdo_j9, 



bb' 



(E.l) 
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X (dbda^RdJ ■ ■ ■ {da^Rd^)9bb' , 

{dbRd){didi,Rd){da,RdJ ■ ■ ■ {da^Rd^){,db'Rdo-J 



(E.2) 



^d. 



■D 



^di-dD^ai-aj^ 



X {9bb9b'b' + 9bb'9bb' + 9bb'9bb'), 



(E.3) 



which we consider in turn. 

The first term contains the ingredient 



(E.4) 



which is identical to the term in Eq. ([HI), computed in Section [221 By substituting its 
value into Eq. (1E.1|) . we obtain 



It is straightforward to observe (via partial integration) that J'^ is identical to J'^, up 
to a sign. As for J'j^ , it is proportional to 



interchange of d and rfj. With the only possible values for both d and di being in the 
set {1, . . . , Ip}, J'^ vanishes. 

Appendix F. Velocity for a weak distortion of a hyper-planar vortex 

In this appendix, we consider weak distortions of a ^-dimensional hyper-planar vortex, 
and determine the velocity that the corresponding superflow confers on such distortions. 
We begin with Eq. f[52|) . which is Eq. pUj) but with the replacement Eq. f[5T|) . which 
yields the velocity f/(cr) of points a on the vortex in terms of H. We then expand 
U in powers of H, which we take to be small, and, for the sake of compactness, we 
temporarily omit an overall factor of 2'kIJ)T (-0/2) /Ati^/"^ which multiplies the velocity. 
Hence, we arrive at the formula 



Jl^ = Ip\ ^/detg Ndj,_,dD {dbdb'Rdo-i) 9bb' ■ 



(E.5) 





(F.l) 




X 



in^) - m)d 

\F{a)-F{s)\ 



d.D-1 



D 
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{H{a) - H{s)),^_^ 



\F(a) - FisM"" 



+ ed,-dD j d^sda,Fd,is) ■ --da^Fd^is) 
- Dedi-do J d^sda^Fd^is) ■ --daj^Fd^is) 

{F{a)-F{s)). 

+ O [H^] . (F.2) 

As F{s) describes a hyper-plane, its derivatives with respect to Sa are simply Kronecker 
deltas: daFd{s) = 6ad- Furthermore, the term of zeroth order in H on the right hand 
side of Eq. (Q, viz.. 



Gdi- do J Sa^^di ■ ■ ■ Saj^c 



\F(a) - FisM"" 



cr - s 



J \F{a)-F{s)\ 

vanishes because a^, can only take values between ai and and the term is already 
skew-symmetric in these indices. The last term of first order in H on the right hand 
side of Eq. (IF. 20 also vanishes, via the orthogonality of F and H. 
We observe that the remaining terms for the velocity, i.e.. 



„ 



(o- - s)„, 5, 



u'do- 



um'=1 



\F(a)-F(s)f 



+ Jd s e,,...n-2d,_,d, \F{a)-Fisr ' ^ 
are nonzero only for the components t/a-i and Ur,. In order to keep the notation 
compact, it is useful to introduce a new index, 7, which takes only the values D — 1 and 
D. As a last step, we rearrange the indices of the Levi-Civita symbols, arriving at the 
{D — 1)*'^ and D^^ components of the velocity, to first order in the height: 

d^s 



7=D-1 



((if (a) - His))^ - J2(duH,is))ia - s)J) , (F.4) 



X 

with all other components of U vanishing to this order. 

Appendix G. Fourier transform of the velocity for the weak distortion of a 
hyper-planar vortex 



The required Fourier transform is of the convolution on the right-hand side of Eq. 
and is of the form 

e-- / d'^s (--'^^(--'^^ ph. (G.l) 
J \(y-n ds^dsp 
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(G.3) 



(G.2) 



The remaining integral is readily computed in hyper-spherical coordinates: 



J d^ae-''^-" 




(G.4) 





/ 

^0 



■TT 



,0-1 



d6 sin 



0-2 



9 {s cos eye 



2-h 



iqs cos 9 



(G.5) 



r(W2) 



in(i/|g|0, 



(G.6) 



where we have imposed a short- distance cutoff at s = ^ to render the integral convergent. 
The spatial Fourier transform evaluated in the present appendix, together with the 
temporal one, convert the differential (in time), integro-differential (in space) equation 
of motion fl55|) into algebraic form; see Eq. fl59|) . 
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